Second wind of the Dulong-Petit Law at a quantum critical point 
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Renewed interest in 3 He physics has been stimulated by experimental observation of non-Fermi- 
liquid behavior of dense 3 He films at low temperatures. Abnormal behavior of the specific heat 
C(T) of two-dimensional liquid 3 He is demonstrated in the occurrence of a T-independent /3 term in 
C(T). To uncover the origin of this phenomenon, we have considered the group velocity of transverse 
zero sound propagating in a strongly correlated Fermi liquid. For the first time, it is shown that if 
two-dimensional liquid 3 He is located in the vicinity of the quantum critical point associated with 
a divergent quasiparticle effective mass, the group velocity depends strongly on temperature and 
vanishes as T is lowered toward zero. The predicted vigorous dependence of the group velocity can 
be detected in experimental measurements on liquid 3 He films. We have demonstrated that the 
contribution to the specific heat coming from the boson part of the free energy due to the transverse 
zero-sound mode follows the Dulong-Petit Law. In the case of two-dimensional liquid 3 He, the 
specific heat becomes independent of temperature at some characteristic temperature of a few mK. 
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Almost two hundred years ago, Pierre-Louis Dulong 
and Alexis-Therese Petit [l[ discovered experimentally 
that the specific heat C(T) of a crystal is close to con- 
stant independent of the temperature T. This behavior, 
attributed to lattice vibrations — i.e. phonons — is known 
as the Dulong-Petit Law. Later, Ludwig Boltzmann [2} 
reproduced the results of Dulong and Petit quantitatively 
in terms of the equipartition principle. However, subse- 
quent measurements at low temperatures demonstrated 
that C(T) drops rapidly as the T is lowered toward zero, 
in sharp contrast to Boltzmann's theory. In 1912, Pe- 
ter Debye [3] developed a quantum theory for evaluation 
of the phonon part of the specific heat of solids, cor- 
rectly explaining the empirical behavior C(T) ~ T 3 of 
the lattice component as T — > 0. In the Debye theory, 
the T-independence of C(T) is recovered at T > To, 
where Tjj is a critical temperature corresponding to the 
saturation of the phonon spectrum. With the advent of 
the Landau theory of quantum liquids [4], predicting a 
linear- in-T dependence of C (T) for the specific heat con- 
tributed by itinerant fermions, our understanding of the 
low-temperature thermodynamic properties of solids and 
liquids thus seemed to be firmly established. However, 
recent measurements of the specific heat of two- 

dimensional (2D) 3 He as realized 3 He films absorbed on 
graphite preplated with a 4 He bilayer, reveal behavior 
strongly antithetical to established wisdom, which calls 
for a new understanding of the low-temperature thermo- 
dynamics of strongly correlated many-fermion systems. 

Owing to its status as a fundamental exemplar of the 
class of strongly interacting many-fermion systems, liquid 
3 He remains a valuable touchstone for low-temperature 
condensed-matter physics. In recent years, interest in 
3 He physics has been driven by the the observation of 
non-Fermi-liquid (NFL) behavior of dense 3 He films at 



the lowest temperatures T ~ 1 mK reached experimen- 
tally [Bl-fll"j|. In particular, measurements of the specific 
heat C(T) in the 2D 3 He system show the presence of 
a term (3 tending to a finite value as T — > 0. Such 
behavior contrasts sharply with that of its counterpart, 
three-dimensional (3D) liquid 3 He; for this system, the 
lower the temperature, the better Landau Fermi-liquid 
(FL) theory works. (Here we shall not consider super- 
fluid phases of 3 He.) 




FIG. 1: Phase diagram of the 2D liquid 3 He system. The 
region defined by z — p/pfc < 1 is divided into LFL and NFL 
domains separated by a solid line. The dependence M*(z) oc 
(1 — z)" 1 is shown by the solid line approaching the dashed 
asymptote, thus depicting the divergence of the effective mass 
at the quantum critical point (« = 1, T = 0) indicated by the 
arrow. In the region z > 1, the fermion condensate (FC) sets 
in and Dulong-Petit behavior of the specific heat is realized 
for the strongly correlated quantum many-fermion system (as 
represented by the dash horizontal line at T — 0). 

In seeking the origin of the anomalous contribution (3 
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remaining in C(T) at the lowest temperatures attained, 
it is instructive to examine the schematic low-T phase 
diagram of 2D liquid 3 He shown in Fig. 1. The essential 
features of this picture are documented by the cited ex- 
periments on 3 He films. The effective coupling parameter 
is represented by z — p/poo, where p is the number den- 
sity of the system and p^ is the critical density at which 
a quantum critical point (QCP) occurs. This QCP is as- 
sociated with a divergence of the effective mass M*(z), 
portrayed in Fig. 1 by the curve (in red on line) that 
approaches the dashed asymptote at z = 1. The T — z 
phase plane is divided into regions of 2D Fermi-Liquid 
(FL) and non-Fermi-Liquid (NFL) behavior. The part of 
the diagram where z < 1 consists of a FL region at lower 
T and a NFL region at higher T, separated by a solid 
curve. The regime where z > 1 belongs to a NFL state 
with specific heat taking a finite value fi(p) at very low 
temperatures. The physical source of this excess heat ca- 
pacity has not been established with certainty, although 
it is supposed that the j3 anomaly is related to peculiar- 
ities of the substrate on which the 3 He film is placed. 

In this letter we propose that the observed /3 term in 
C(T) can instead have its origin in an intrinsic mecha- 
nism analogous to that producing the classical Dulong- 
Petit behavior in solids. It is shown for the first time that 
in systems (such as 2D 3 He) containing a fermion con- 
densate (FC), the group velocity of transverse zero sound 
depends strongly on temperature. It is this dependence 
that gives rise to the (3 term, granting the Dulong-Petit 
Law a "second wind." 

As indicated above, the most challenging feature of 
the NFL behavior of liquid 3 He films involves the specific 
heat C(T). According to Landau theory, C(T) varies lin- 
early with T, and at low film densities the experimental 
behavior of the specific heat of 2D liquid 3 He is in agree- 
ment with FL theory. However, for relatively dense 3 He 
films, this agreement is found to hold only at sufficiently 
high temperatures. If T is lowered into the millikelvin 
region, the function C(T) ceases to fall toward zero and 
becomes flat [ISED]. 

The common explanation 0, 0, [l2[ of the flattening of 
C(T) seen in these experiments imputes the phenomenon 
to disorder associated with the substrate that supports 
the 3 He film. More specifically, it is considered that 
there exists weak heterogeneity of the substrate (steps 
and edges on its surface) such that quasiparticles, be- 
ing delocalized from it, give rise to the low-temperature 
feature f3 of the heat capacity [6j. Even if we disregard 
certain unjustified assumptions made in Ref. H3, there 
remains the disparate fact that the emergent constant 
term in C(T) is of comparable order for different sub- 
strates JE S El- Furthermore, the explanation posed 
in Ref. [12| implies that the departure of C(T) from FL 
predictions shrinks as the film density increases, since ef- 
fects of disorder are most pronounced in weakly interact- 
ing systems. Contrariwise, the anomaly in C(T) makes 



its appearance in the density region where the effective 
mass M* is greatly enhanced [y, \ll\ and the 2D liquid 
3 He system becomes strongly correlated. This reasoning 
compels us to consider that the NFL behavior of C(T) 
is an intrinsic feature of 2D liquid 3 He, which is asso- 
ciated with the divergence of the effective mass rather 
than with disorder. The flattening of the curve C(T) as 
seen in 3 Hc films is by no means a unique phenomenon. 
Indeed, as expressed in the Dulong-Petit (DP) law, the 
specific heat C(T) of solids remains independent of T as 
long as T exceeds the Debye temperature fiu, which is 
determined by the saturation of the phonon spectrum of 
the crystal lattice. Normally, the value of SId is suffi- 
ciently high that the DP law belongs to classical physics. 
However, we will argue that the DP behavior of C(T) 
can also make its appearance at extremely low temper- 
atures in strongly correlated Fermi systems, with zero 
sound playing the role of phonons. 

To clarify the details of this phenomenon and calculate 
the specific heat C (T) , we evaluate a part Fb of the free 
energy F associated with the collective spectrum uj(k) = 
ck, based on the standard formula 
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J^-^—Jlm[lnD-\k,u)]d V , (I) 



where D(k,ui) is the boson propagator, and dv is an el- 
ement of momentum space. Upon integration by parts 
this formula is recast to 



f duj , / , /T \ f ( dD- 1 (k,io)ldw\ , 



(2) 

If damping of the collective branch is negligible (the 
case addressed here), then D^ 1 (k,u) ~ (u — ck) and 
dD- 1 (k,u))/du ~ 1, while ImD^ik,^) ~ 5(u - ck), 
and we arrive at the textbook formula 



F B =T J In (l - e- cfe / T ) 9(Q ~ ck)di 



(3) 



where Qq is the characteristic frequency of zero sound. 
At T > n , the factor ln(l - e~"/ T ) reduces to ln(w/T), 
yielding the result 



F B (T)<xT\n(Slo/T), 



(4) 



which, upon the double differentiation, leads to the DP 
law C(T) = const. At first sight, this law has nothing to 
do with the situation in 2D liquid 3 He. Its Fermi energy 
e a F is around I K at densities where the Sommcrfcld ratio 
C(T)/T soars upward as T — > 0, while fio must be lower 
than T ~ 1 mK. Indeed, in any conventional Fermi liquid, 
including 3D liquid 3 He, there is no collective degree of 
freedom whose spectrum is saturated at such low ratios 
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This conclusion remains valid assuming 2D liquid 3 He 
is an ordinary Fermi liquid. However, as seen from Fig.[T] 
if the quantum critical point is reached at T — > and 
some critical density where the effective mass M* (p) 



diverges, as it does in the present case [6H8|, 111], |l3j, the 
situation changes dramatically. This is demonstrated 
explicitly in the results of standard FL calculations of 
the velocity Ct of transverse zero sound, which satisfies 

mm 



Ct , c t + V F 
- — In 

2v F c t - v F 



- 1 = 



Fi-6 



3F 1 (c 2 /v 2 



(5) 



where vp = pp/M* is the Fermi velocity and F± = 
PfM* fi/ir 2 is a dimensionless version of the Landau first 
harmonic /i [3, H E3] • The divergence of the effective 
mass M* at the QCP implies that at the critical density 
determined by fippM/n 2 = 3 la. 17|. one has 



! (P) 



P F 



5M*(p)M 



0. 



(6) 



whereas the sound velocity c s remains finite in this limit 

We see then that in case the effective mass M* di- 
verges, the group velocity c t vanishes as y/l/M*. Flat- 
tening of the single-particle spectrum e(p) prevails as 
long as \p — pf\/pf < M/M*, implying that the trans- 
verse mode softens only for rather small wave numbers 
k ~ pfM/M*. Unfortunately, the associated numerical 
prefactor cannot be established, rendering estimation of 
~ (pp/M)-^/M/M* correspondingly uncertain. Nev- 
ertheless, one cannot exclude a significant enhancement 
of the Sommerfeld ratio C(T) /T at T ~ 1 mK due to 
softening of the transverse zero sound in the precritical 
density region. 

At T — > and densities exceeding p^, the system 
undergoes a cascade of topological phase transitions in 
which the Fermi surface acquires additional sheets 
211 ] . As indicated in Fig. [U FL theory continues to hold 
with quasiparticle momentum distribution n(p) satisfy- 
ing n 2 = n, until a greater critical density ppc is reached 
where a new phase transition, known as fermion conden- 
sation, takes place (2l| - [27| . Beyond the point of fermion 
condensation, the single-particle spectrum e(p) acquires 
a flat portion. The range L of momentum space adjacent 
to the Fermi surface where the FC resides depends on the 
difference between the effective coupling constant and its 
critical value. As will be seen, L is a new dimensional 
parameter that serves to determine the key quantity ^Iq. 

At finite T, the dispersion of the FC spectrum e(p) 
existing at p > P fc acq uires a nonzero value proportional 
to temperature [21, 24- 2a): 



e(p,T) = T\n 1 ^\ Pl <p< Pf , (7) 

where < < 1 is the FC momentum distribution 

and pi and p / are the lower and upper boundaries of the 



FC domain in momentum space. Consequently, in the 
whole FC region, the FC group velocity, given by 



v(p,T) 



de(p) 
dp 



<9n* (p) I dp 
n*(p)(l - n*(p)) ' 



= -T — - ; 7." " —, — , p t <p<p f , 

(8) 

is proportional to T. Significantly, in the density interval 
Poo < P < Pfc the formula describes correctly the 
single-particle spectrum e(p, T) in case the temperature 



T exceeds a very low transition temperature [21[ . The FC 



itself contributes a T-independent term to the entropy 
S; hence its contribution to the specific heat C(T) — 
TdS/dT is zero. Accordingly, we focus on the zero-sound 
contribution to C(T) in systems having a FC. 

Due to the fundamental difference between the FC 
single-particle spectrum and that of the remainder of 
the Fermi liquid, a system having FC is, in fact, a two- 
component system. Remarkably, the FC subsystem pos- 
sesses its own set of zero-sound modes, whose wave num- 
bers are relatively small, not exceeding L — (pf—Pi) > 0. 
The mode of prime interest for our analysis is that of 
transverse zero sound. As may be seen by comparison of 
formulas ([5]) and ([5)1, its velocity Ct depends on tempera- 
ture so as to vanish like y/T as T — » 0. 

To verify the latter property explicitly, we observe first 
that for systems with a rather small proportion of FC, 
evaluation of the spectrum of collective excitations may 
be performed by employing the familiar FL kinetic equa- 
tion 0, m 



(w — kv) 5n(p) = — kn 



dn(p) 
dp 



J r (p,Pi)5n(pi)dwi. 



(9) 

Focusing on transverse zero sound in 2D liquid 3 He, one 
need only retain the term in the Landau amplitude F 
proportional to the first harmonic f\. To proceed fur- 
ther we make the usual identification (c< — cos8)Sn(p) = 
(dn(p) I dp) 4>(n) , where cos# = kv/kv. Equation (|9]) 
then becomes 



-fxp F cos I cosx 



dn(pi)/dpi 
c t - v(T) cos( 



dpidOi 



" (2tt)2 ' 
(10) 

where cos \ — cos cos 0\ +sin 9 sin Q\ , while v(T) is given 
by equation ([SJ. The solution describing transverse zero 
sound is </>(n) ~ sin 9 cos 9. 

We see immediately that Ct v(T) ~ T; therefore 
the transverse sound in question does not suffer Landau 
damping. In this situation, we are led to the simple result 



2 p F f dn(p) 

c * = -smJ -dp~ v ^ dp ^ 

upon keeping just the leading relevant term v(T) cos9/c 2 
of the expansion of (ct — f(T) cos and executing 
straightforward manipulations. Factoring out an aver- 
age value of the group velocity v(p, T) oc T/pp, we arrive 
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at the stated behavior 



c t (k) 




(12) 



for wave numbers k not exceeding the FC range L. Trans- 
verse sound can of course propagate in the other, noncon- 
densed subsystem of 2D liquid 3 He consisting of quasi- 
particles with normal dispersion 14, lj| 28 1. However, 
its group velocity is T-independent, so the corresponding 
contribution to the free energy is irrelevant. 

As noted above, the characteristic wave number of the 
soft transverse zero-sound mode is given by the FC range 
L{p) = pf — Pi, treated here as an input parameter. The 
key quantity Qq is therefore estimated as fio — k max ct, 
where fc max is the maximum value of the zero sound mo- 
mentum at which the zero sound still exists. In our case 
the zero sound is associated directly with the FC; hence 
fcmax — VLpf and we have 



r» 



TLp f 
M 



(13) 



As long as the inequality Lpp/M < T is met (or equiva- 
lently, T/e° F > L/p F holds), the ratio Qq/T is small, and 
we are led to the DP result C(T) = const. Then, in spite 
of the low temperature, C behaves as if the system were 
situated in the classical limit rather than at the QCP. 
Such a behavior is ensured by the fact that the system 
contains a macroscopic subsystem with heavy quasiparti- 
cles. As the temperature ultimately goes down to zero at 
the fixed density p, the inequality Lpp/M < T eventually 
fails, the quantum regime is restored and the dominant 
contribution to C comes from the "normal" fermions. In 
other words, there exists an extremely low temperature 
To below which the usual FL behavior of zero sound is 
recovered. 

Interestingly, the value of the constant term in C(T) 
can be evaluated in closed form in terms of the FC range 
L. Upon inserting Wt(fc) = c t k into Eq. (j3|) and integrat- 
ing, the T-independent term in the specific heat is found 
to be 



C_ 

N 



Lpf 
8irp 



(14) 



where N is the number of atoms in the film. The FC 
range parameter L also enters the result derived analo- 
gously for the spin susceptibility \- The FC component 
of x is given by 21. 29| 



X.(F) ~ xc(T)—, 

PF 



(15) 



where xc(T) = [i 2 B p/T. 

The results (|14[) and ([151) jointly establish an un- 
ambiguous relation within our model between the T- 
independent term in the specific heat C (T) and the Curie 



component of the spin susceptibility x{T) (which has also 
been observed experimentally @,H[)- This relation can 
be tested using existing experimental data @. The T- 
independent specific heat C/N exists in the density re- 
gion around p = 9.5 nm~ 2 . Being referred to one parti- 
cle, it is readily evaluated from (3 ~ 0.25 mJ/K. One finds 
C/N ~ 0.01, yielding L/p F ^ 0.05. On the other hand, 
the data for the spin susceptibility given in Fig. 2(B) 
of Ref. @ supports a Curie-like component at p = 9.25 
nm~ 2 . The value of the corresponding numerical factor 
extracted from the data, which according to Eq. (Ti"5"j) is to 
be identified with the ratio L/pp, is approximately 0.07. 
Given the uncertainties involved, our model is consistent 
with the experimental data reported in [6]. 

In summary, we have analyzed the group velocity of 
transverse zero sound propagating in a strongly corre- 
lated Fermi liquid. We have shown for the first time that 
if two-dimensional liquid 3 He is located in the vicinity 
of the quantum critical point associated with a diver- 
gent quasiparticle effective mass, the group velocity de- 
pends strongly on temperature and vanishes at diminish- 
ing temperatures. Such a vivid dependence of the group 
velocity can be detected in experimental measurements 
on the liquid. We have demonstrated that the contri- 
bution to the specific heat coming from the boson part 
of the free energy contributed by the transverse zero- 
sound mode follows the Dulong-Petit law. Accordingly, 
the specific heat becomes independent of temperature 
at some characteristic temperature. In the case of two- 
dimensional liquid 3 He, this temperature can be a few 
mK. At sufficiently lower temperature the usual FL be- 
havior of zero sound is recovered. The model developed 
from the analysis is found to be in reasonable agreement 
with experimental measurements. 
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